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ABSTRACT 

Recent observations of the white dwarf (WD) populations in the Galactic globular cluster NGC 6397 
suggest that WDs receive a kick of a few knis^^ shortly before they are born. Using our Monte Carlo 
cluster evolution code, which includes accurate treatments of all relevant physical processes operating 
in globular clusters, we study the effects of the kicks on their host cluster and on the WD population 
itself. We find that in clusters whose velocity dispersion is comparable to the kick speed, WD kicks 
are a significant energy source for the cluster, prolonging the initial cluster core contraction phase 
significantly so that at late times the cluster core to half-mass radius ratio is a factor of up to ^ 10 
larger than in the no-kick case. WD kicks thus represent a possible resolution of the large discrepancy 
between observed and theoretically predicted values of this key structural parameter. Our modeling 
also reproduces the observed trend for younger WDs to be more extended in their radial distribution 
in the cluster than older WDs. 

Subject headings: globular clusters: general — methods: numerical — stellar dynamics — white dwarfs 



L WHITE DWARF BIRTH KICKS 

Recent i/ST observations of the Galactic globular clus- 
ter NGC 6397 reveal that its youngest white dwarfs 
(WDs) are significantly more extended in their radial 
distribution within the cluster t han their older visible 
counterparts (jDavis "eraI1l2008br i. Since the progenitors 
of the WDs are some of the most massive stars in the 
cluster at ^ 0.8 Mq, and since the older WDs (of mass 
~ 0.5 Mq) have had time to diffuse to larger radii than 
where they formed due to mass segregation, one would 
naively expect the opposite — namely that the youngest 
WDs should be more centrally concentrated than their 
older counterparts. One possible interpretation is that 
WDs rece ive a kick of a few kms~^ shor t ly befo re they 
are born (|Davis et all l2008bt [Hevl [20071 . 12008bl lah. In- 
deed, in open star clusters, where the velocity dispersion 
is ~ 1 kms~^ and the putative WD kick would be enough 
to eject it from the cluster, there is an observed relative 
dearth of WDs (iFellhauer et al.l l2003l: iKahrai et alll2001l : 
IWeidemannlll977f K 

The mechanism generating WD "natal" kicks is likely 
asymmetric mass loss late in the asymptotic giant branch 
(AGB) phase. Although AGB winds have not been di- 
rectly observed, the effects of the resulting mass loss have 
been, and there a re good theoretical reason s to expect 
such winds exist (jVassiliadis fc WoodI |1993[ ). Alterna- 
tively, the kick could be produced result of small 
asymmetry during the helium core flash (Ivanova. pri- 
vate communication). Independent of precisely when in 
the evolution of a star the kick occurs, the observed rota- 
tion rates of WDs arc consistent with non- axisyrnmetri c 
mass loss at some point in their evolution ()Spruitj|1998D . 
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Analogous to the hydrogen-burning main sequence in 
stars, star clusters are thought to eventually enter a long 
lived binary-burning phase in which inelastic dynamical 
scattering interactions of binaries with other stars gener- 
ate enough e nergy to prevent the clust er core from col- 
lapsing (e.g.. iBinnev fc Tremaind [20081 ). The physics of 
the binary-burning phase has been studied intensively. It 
is only recently, however, that different numerical meth- 
ods have begun to agree on the structural properties 
of a cluster in the binary burning phase (at least in 
the e qual-mass case; iHeggie et aLll2006HRe geau fc Rasiol 
l2007t l. Unfortunately, the predicted value of the core 
to half-mass radius ratio rc/vh in the binary-burning 
phase is at least ~ 10 times smaller than what is ob- 
served for the majority of Galactic globular clusters (e.g.. 
iFregeau fc Rasiol |2007[ ). Although differences in the ob- 
servational and theoretical definitions of the core radiu s 
may account for part of the difference (jHurlevI l2007t l , 
there is still a significant discrepancy. If the bulk of 
Galactic clusters are indeed currently in the binary- 
burning phase, a core radius larger than expected from 
binary burning suggests the existence of an alternate 
core energy source. In this vein, it has been suggested 
that perhaps tens of Galactic clusters contain a central 
intermediate-mass black hole (IMBH; Trenti i2006l) , pow- 
ering their cores to the observed sizes. Or the required 
energy could come from expedited stellar evolution due 
to stellar collisions (|Chatteriee et al.l l 2007D. or prolonged 
mass segregation of compact objects (jMerritt et al.l[20M : 
iMackev et al.|[2007t ) . Alternatively, it could be that most 
clusters are simply not currently in the binary burning 
phase, in which case their init ial properties may play an 
important role (lFregeaull2008f). 

Since for a standard iKroupa et al.l (jl993f ) initial mass 
function, roughly 15% of stars in a cluster will become 
WDs within 12 Gyr, a WD birth kick of a few kms~^ 
represents a potentially significant core energy source 
(roughly 15% percent of the total cluster energy for a 
WD kick speed comparable to the velocity dispersion) — 
one that could resolve the discrepancy between theory 
and observations on the structural parameters of clus- 
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ters. In this paper wc consider the dynamical effects of 
WD birth kicks by simulating the evolution of star clus- 
ters in which WDs receive kicks of a few kms~^ at birth. 
Our simulation method treats all relevant physics in clus- 
ters self-consistently, as described below. Wc find that 
WD birth kicks can significantly prolong the initial clus- 
ter evolutionary phase of core contraction, resulting in 
core sizes at late times that are larger than for clusters 
without kicks. We also find that birth kicks produce an 
observational signature that is consistent with the ob- 
served radial distributions of WDs in NGC 6397. 

2. SIMULATIONS 

To simulate the evolution of clusters with WDs 
that experience birth kicks, we use the well-tested 

Henon M onte Carlo metho d, as described in detail 

elsewhere (iJoshi et all l2000L l200H iFregeau et all 120031 : 
iFregeau fc Rasioll2007l [20091) . Briefly, it treats the effects 
of two-body relaxation in the orbit-averaged approxima- 
tion and, since it provides for a physical realization of 
the cluster at each timestep, has been extended to in- 
clude a treatment of all other relevant physics, including 
strong dynamical interactions of binaries, physical stel- 
lar collisions, stellar evolution of single stars and binaries, 
and the effects of a Galactic tide. It has shown excellent 
agreement with direct iV-body methods for the evolution 
of clus ters with dynamically si gnificant populations of bi- 
naries (|Fregeau fc Rasiall2007D , and tho se in which the ef- 
fects of stellar evolution are dominant (|Fregeau fc Rasiol 
l2009t ) . To model the effects of WD birth kicks we simply 
give a randomly oriented velocity kick to a star as soon 
as it becomes a WD. For simplicity, the kick is set to a 
fixed value between 2 and 9 kms~^ for each simulation. 

Modeling individual clusters wit h any degree of de- 
tail is a difficult endeavor (e.g., iHurlev et al.l l2008l : 
iGiersz fcll eggie 2009). While we do not attempt to 
model NGC 6397 in this paper, we find that our mod- 
els do have roughly the same core velocity dispersion at 
late times, and the model with 4kms~^ WD kicks has 
roughly the same ratio of core to half-mass radius at 
very late times (to within a factor of 2 of the value in 
iNovola fc Gebhardtll2006[ ). Our standard initial model 
is a Wo = 7.5 King model with a virial radius of 5pc, 
consisting of = 3 x 10^ objects, 1% of which are 
binaries (similar to NGC 6397's inferred initial binary 
fraction; iDavis et ~al.. 2008a ) . No initial mass segregation 
is assume d. Single stars and binary primaries are cho- 
sen from a lKroupa et al.l (|l993f ) distribution from 0.15 to 
18.5 M0. Binary secondaries are chosen from a flat mass 
ratio distribution with a minimum at 0.15 A/q, semima- 
jor axes are chosen uniformly in the log from 5(i?i -t- R2) 
(where Ri are the stellar radii) to the hard-soft bound- 
ary, and eccentricities are chosen from a thermal distri- 
bution truncated at the upper end to prevent pericenter 
approaches smaller than 5(R^ + R2)- Clus t er me tallic- 
ity is set to Z = 0.001, as in iHurlev erall ()2008f ). Al- 
though the metallicity of NGC 6397 is a fa ctor of ~ 5 
smaller, as described in IHurlev et al.l (j2008f ) stellar evo- 
lutionary processes don't depend strongly on metallicity 
below Z « 0.001. The Galactic tide represents that of a 
point-mass host galaxy of mass 10^^ Mq, with the cluster 
in a circular orbit at 6.9 kpc. 

Figure [Tj shows the evolution of Vc/rh for our standard 
model with no WD birth kicks, 4kms^^ birth kicks, and 



6kms~^ birth kicks. Around llGyr, the model with- 
out kicks enters the binary burning phase, during which 
time rc/rh « 0.01. (Note that we use the standard 
three-dimensional, A^-b ody definition of the core radius; 
IFregeau fc Rasiol l2009f ). The models with kicks show no 
evidence of reaching a binary burning phase in their 16 
Gyr of evolution. Instead, they show a prolonged, grad- 
ual contraction of the core with time, si milar to the be- 
havior found in A^-body simulations by iHurlevI ([2007) . 
As expected, the model with 6kms~^ birth kicks dis- 
plays an appreciably larger core radius over time com- 
pared to the 4kms~^ kick model. At late times, the 
models with kicks display large values of rc/rh, with val- 
ues up to ~ 10 times larger than the model without kicks. 
We note that for the model with 4kms~^ kicks we find 
a value of Vc/rh k, 0.11 at an age of 12 Gyr, which is sig- 
nificantly large r than the observed value of Vc/rh ~ 0.03 
for NGC 6397 (iNovola fc Gebhardtll200l . The 4kms-i 
kick model reaches values of Tc/rh ~ 0.05 at unphysi- 
cally late times, but further simulations are required to 
determine if a better match to NGC 6397 can be made 
with a suitable choice of initial conditions. 

Figure[2]shows the cumulative three-dimensional radial 
distributions of WDs at 12 Gyr in the models with no 
kick and with 6kms^^ kicks, broken down into "young" 
and "old" populations. The young population consists 
of those WDs with L > 1.2 x 10~"^ L p, (correspond- 
ing t o ages < 4 Gyr for a 0.5 Af© WD; iBergeron et al.l 
Il995f ). while the old population consists of WDs with 
1.9 X 10~^ < L < 7.9 X 1O~^L0 (correspon ding to ages 
between 6 and 10 Gyr for a 0.5 Mq WD). In lDavis et al.l 
(|2008bf) analogous populations are chosen with the re- 
striction that the young WDs have ages less than ^ 3 
times the local relaxation time of 0.29 Gyr, so that mass 
segregation has had minimal effect on their radial dis- 
tribution. The e ffects of mass segre gation decay expo- 
nentially in time (|Fregeau et al.ll2003 ) , so they should be 
apparent for any division between young and old WDs. 
To maximize the evidence of the effect, we take the cut 
at roughly a half-mass relaxation time, which is w 5 Gyr 
in our models at an age of 12 Gyr. The top panel shows 
the cumulative radial distribution of young and old WDs 
in the model with no WD birth kicks. The young pop- 
ulation is statistically significantly more centrally con- 
centrated than the old population, as might be expected 
when one considers that the progenitors of the young 
WDs are ~ 0.8 Af© stars and hence more massive than 
the ^ 0.5 A^o "old" WDs which have had time to mass 
segregate out to their larger cluster radii. The bottom 
panel shows the cumulative distributions for the model 
with 6 kms^^ kicks. The young population is clearly and 
significantly more radially extended than t he old popu- 
lation. This is consistent with the results of iHevll (|2007l . 
2008bh. and of course with the observational results of 
Davis et al.l 'l2008bD. 

iDavis et al.l (|2008b[ ) observed the positions of the WDs 
in their sample in projected radius between ~ 3 and 
~ 6pc (1.7 and 3.5rft). Figure [31 shows the cumula- 
tive projected radial distributions at ~ 12 Gyr between 
12 and 24 pc (1.7 and 3.5r/i in our model) of old and 
young WDs in the model with 4kms~^ WD birth kicks. 
The young WDs are clearly more radially ex tended than 
the old population, in agreement with what iDavis et al.l 
(|2008bf l find in NGC 6397. Of course, while observers 
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Fig. 1. — Evolution of core to half-mass radius ratio, rc/r/i, for 
models with no WD birth kicks, 4kms^^ birth kicks, and 6kms~^ 
birth kicks. 

can only take snapshots of stellar positions at the cur- 
rent time, with numerical modeling we have the lux- 
ury of making many statistical realizations of a model. 
Figure [3] shows one of many possible realizations — some 
yield young WD populations that are even more radi- 
ally extended, some less. In no case where we looked 
did we find one where the old population is statistically 
significantly more radially extended than the young pop- 
ulation. We note that the average mass of the young WD 
population plotted in this figure is 0.55 and that of 
the old population is 0.56 Mq, confirming that mass seg- 
regation cannot explain the observed difference in radial 
distributions. For Figure [5] the average masses in the 
young and old WD populations are 0.558 and 0.563 M©, 
respectively, for the top panel and 0.55 and 0.57 Mq for 
the bottom panel. 

In attempting to make quantitative conclusions about 
the WD kick speed, we are faced with the lack of detailed 
agreement between our model and NGC 6397. In the 
cursory mapping of parameter space we've done for this 
paper (consisting of ~ 50 simulations), we've found that 
even very small kicks of 2 kms^^ are sufficient to "puff 
up" the core noticeably. On the other hand, the pro- 
jected radial distribution of WDs in the same field as the 
observations and with the same luminosity cuts is some- 
what more sensitive to the kick speed — too low and the 
effect is not apparent, too high and too few WDs remain 
for meaningful statistical comparisons. Without getting 
mired in detailed modeling we simply note that, as seen 
in Fig. [21 an extended radial distribution for the young 
WDs compared to the older WDs is a robust consequence 
of WD birth kicks. Furthermore, a kick spe ed of '-^ 3 to 
5 kms~ ^, as inferred from observations bv iDavis et al.l 
(|2008br ). is sufficient to prolong the initial evolutionary 
process of core contraction and result in a core radius 
at the current time that is significantly larger than that 
expected from simple binary burning. 

3. DISCUSSION 

Through realistic numerical simulations of globular 
cluster evolution, we have shown that WD birth kicks 
of a few kms~^ act as an energy source that can "puff 
up" a cluster's core potentially significantly. In a cluster 
whose velocity dispersion is comparable to the WD kick 
speed, Tc/r/i can be increased by a factor of up to 10 
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Fig. 2. — Cumulative three-dimensional radial distributions of 
old (solid curve) and young (dashed curve) WD populations at 
~ 12 Gyr in models without WD birth kicks (top panel), and with 
kicks of 6 km (bottom panel). The young population consists of 
those WDs with L > 1.2 X 10~* Lq (corresponding to ages < 4 Gyr 
for a 0.5 Mq WD), while the old population consists of WDs with 
1.9 X 10~^ ^ L < 7.9 X 10~^Lq (corresponding to ages between 6 
and 10 Gyr for a 0.5 Mq WD). The probability for obtaining the 
two-sided K-S statistic for each pair of distributions is also shown. 
In the no kick case the young population is statistically significantly 
more centrally concentrated than the old population, while in the 
6kms~^ kick case it has a more extended radial distribution. 

at late times relative to the no-kick case. Coincidentally, 
theoretical predictions of Vc/rh for a cluster whose core 
is supported against collapse purely by inelastic scat- 
tering interactions of binaries (binary burning) have re- 
cently been shown to be up to a factor of ~ 10 smaller 
than observed for the bulk of Galactic globulars (e.g., 
iFregeau fc Rasioll2007f) . WD kicks offer a natural resolu- 
tion of this large discrepancy for clusters whose velocity 
dispersions are a few kms"-'^, and obviate the need for 
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Fig. 3. — Cumulative projected radial distributions at ~ 12 Gyr 
between 12 and 24 pc (1.7 and S.Srj, in our model) of old (solid 
curve) and young (dashed curve) WD populations in the model 
with 4kms~^ WD birth kicks. WD age cuts are as in the previous 
figure. 



alternate energy sources whose existence remains specu- 
lative, such as IMBHs. 

Although the large core radii of some Galactic globu- 
lars appear t o be naturally explained by central IMBHs 
()Trenti|[2006[ ). there is still only circumstantial evidence 
for their existence in clusters. The evidence that ex- 
ists is indirect (from cluster velocity dispersion pro- 
files), and for M15 and Gl can be explained equally 
well via standard cluster evolution models which con- 
tain a population of neu tron stars or WDs in the core 
(jBaumgardt et al.l l2003al lH). Furthermore, we should 
point out that the numerical models used to infer the 
existence of central IMBHs in,Trcnti ( 2006) contain just 
20,000 particles, and have IMBHs that are a fraction 
M-Qu/Mcius = 0.02 of the cluster mass. Several inde- 
pendent simulations have shown that the mass of the 
stellar progenitor to the IMBH that results from a run- 
away collision scenario (by far the most likely scenario 
for forming an IMBH in a cluster) is generically a frac- 
tion of M-Qu/McXus = 0.002 of the total cluster mass, 
independent of total cluster mass, initial central den- 
sity, properties of the m ass function, and t he degree 



of ini t ial mass segregation (iPortegies Zwart~fc McMillanI 
[200l iGurkan et all [2003 : iFreitag et all l2006bllal) . The 
mass of the IMBH formed is likely to be less than the 
mass of its stellar progenitor, d ue to winds or other mech- 
anisms (jYungelson et al.ll2008[ ). Adopting Mbh/McIus = 
0.002 as an upper hmit, and apply ing the theore t ical sc al- 
ing r^/Th oc (Mbh/McIus)^/'^ from lHeggie et al.l (l2007f) to 
the va lue rdrh ~ 0.3 found in the simulations of iTrentil 
(|200l . the maximum core to half-mass radius ratio for a 
cluster whose core is powered by a central IMBH should 
thus be rc/rh ~ 0.05. This is a ppreciab l y sma ller than 
the value rc/vh ~ 0.3 used by iTrentil (|2006[ ) to infer 
the presence of IMBHs, and is consistent only with the 
~ 20% of clusters with small Vc/rh that are classified 



observationally as core-collapsed. 

We have also shown that the cluster radial distribu- 
tion of WDs with birth k i cks is c onsistent with the ob- 
servations of iDavis et al.l ()2008bf ). While this certainly 
supports the interpretation of kicks as the source of the 
radial distribution, it is far from a direct detection. More 
detailed observations are needed to confirm the existence 
of the kick. In fact, we will soon have relatively precise 
proper motion measurements for the bulk of the WDs 
in the NGC 6397 field, which, with the help of detailed 
simulations, will help pin down the properties of the WD 
kick. Based on our current understanding we can already 
make a few simple predictions, though. There is no rea- 
son to believe the magnitude of the WD kick depends on 
the properties of the host clusters in any significant way. 
We thus expect two measurable consequences. In clus- 
ters whose velocity dispersion is significantly larger than 
the putative kick speed of a few kms~^ (e.g., 47 Tuc), the 
radial distribution of the young WDs should follow or be 
more centrally concentrated than that of the old WDs. 
Observations of WD radial profiles in the high-velocity 
dispersion cluster Omega Cen (~ 20kms~^ in the core ) 
are consistent with this picture (jCalamida et al.ll2008f l. 
In clusters whose velocity dispersion is less than a few 
kms~^ (e.g., NGC 288), there should be a measurable 
paucity of WDs. In fact, in the ~ Ikms^^ velocity dis- 
persion environments o f open clusters, there is thought 
to be a dearth of WDs (|Wilhamsll200l . 

Finally, we note that if the kick does occur during the 
extended asymptotic giant branch (eAGB) phase, there 
are at least two potential complications to our modeling. 
Stellar winds during th e eAGB phase may h ave veloci- 
ties of just ~ lOkms"^ (|Marshall et al.ll2004f ). Since the 
velocity dispersion is comparable to this value (or just 
slightly smaller), there is the very real possibility that 
the stellar wind does not escape the cluster. For clus- 
ters whose velocity dispersion is comparable to the WD 
kick speed, the "heating" of the cluster via the kick is 
direct, and the retained wind mass should not alter our 
results appreciably. However, for clusters whose veloc- 
ity dispersion is less than the kick speed, the "heating" 
is indirect (via cluster mass loss, since the WDs are es- 
caping the cluster), and the retention of the wind mass 
by the cluster may appreciably reduce the "heating" ef- 
fect. Another potential complication is that if the mass 
loss during the eAGB phase takes place on a timescale 
comparable to or longer than the typical stellar orbital 
timescale of ^ 10^ yr, the dynamical result is not that 
of a kick (as we have assumed here), but a weaker, adi- 
abatic modification of the orbit. Mass loss rates of up 
to 10~^MQyr~-'^ have been detected in AGB stars (im- 
plying a mass loss timescale < 10^ yr), but the general 
characteristics of mass loss in this phase are far from 
settled. 



The authors thank B. Hansen and N. Ivanova for help- 
ful discussions. JMF acknowledges support from Chan- 
dra Postdoctoral Fellowship Award PF7-80047. HBR ac- 
knowledges support from The Natural Sciences and En- 
gineering Research Council of Canada, and thanks KITP 
for support during a recent visit during which time this 
paper was conceived. FAR acknowledges support from 
NASA Grant NNG06GI62G and from KITP. This re- 



The Dynamical Effects of White Dwarf Birth Kicks in Globular Star Clusters 



5 



search was completed at KITP and supported in part by the NSF under Grant PHY05-51164. 



REFERENCES 



Baumgardt, H., Hut, P., Makino, J., McMillan, S., & Portegies 

Zwart, S. 2003a, ApJ, 582, L21 
Baumgardt, H., Makino, J., Hut, P., McMillan, S., & Portegies 

Zwart, S. 2003b, ApJ, 589, L25 
Bergeron, P., Wesemael, P., & Beauchamp, A. 1995, PASP, 107, 

1047 

Binney, J. & Tremaine, S. 2008, Galactic Dynamics: Second 
Edition (Princeton University Press) 

Calamida, A., Corsi, C. E., Bono, G., Stetson, P. B., Prada Moroni, 
P. G., Degl'Innocenti, S., Ferraro, I., lannicola, G., Koester, D., 
Pulone, L., Monelli, M., Amico, P., Buonanno, R., Freyhammer, 
L. M., Marchetti, E., Nonino, M., & Romaniello, M. 2008, 
Memorie della Societa Astronomica Italiana (arXiv:0712.0602), 
79, 347 

Chatterjee, S., Fregeau, J. M., & Rasio, F. A. 2007, lAU 

Symposium 246, 2007, Capri, Italy, Poster 
Davis, D. S., Richer, H. B., Anderson, J., Brewer, J., Hurley, J., 

Kalirai, J. S., Rich, R. M., & Stetson, P. B. 2008a, AJ, 135, 2155 
Davis, D. S., Richer, H. B., King, I. R., Anderson, J., Coffey, J., 

Fahlman, G. G., Hurley, J., & KaUrai, J. S. 2008b, MNRAS, 383, 

L20 

Fellhauer, M., Lin, D. N. C, Bolte, M., Aarseth, S. J., & Williams, 

K. A. 2003, ApJ, 595, L53 
Fregeau, J. M. 2008, ApJ, 673, L25 

Fregeau, J. M., Giirkan, M. A., Joshi, K. J., & Rasio, F. A. 2003, 
ApJ, 593, 772 

Fregeau, J. M., Joshi, K. J., Portegies Zwart, S. F., & Rasio, F. A. 

2002, ApJ, 570, 171 
Fregeau, J. M. & Rasio, F. A. 2007, ApJ, 658, 1047 
— . 2009, in preparation 

Freitag, M., Giirkan, M. A., & Rasio, F. A. 2006a, MNRAS, 368, 
141 

Freitag, M., Rasio, F. A., & Baumgardt, H. 2006b, MNRAS, 368, 
121 

Giirkan, M. A., Freitag, M., & Rasio, F. A. 2004, ApJ, 604, 632 



Giersz, M. & Heggie, D. C. 2009, ArXiv e-prints {arXiv:0901.1085) 
Heggie, D. C, Hut, P., Mineshige, S., Makino, J., & Baumgardt, 

H. 2007, PASJ, 59, Lll 
Heggie, D. C, Trenti, M., & Hut, P. 2006, MNRAS, 368, 677 
Heyl, J. 2007, MNRAS, 381, L70 
— . 2008a, MNRAS, 390, 622 
Heyl, J. S. 2008b, MNRAS, 385, 231 
Hurley, J. R. 2007, MNRAS, 379, 93 

Hurley, J. R., Shara, M. M., Richer, H. B., King, I. R., Saul Davis, 
D., Kalirai, J. S., Hansen, B. M. S., Dotter, A., Anderson, J., 
Fahlman, G. G., & Rich, R. M. 2008, AJ, 135, 2129 
Joshi, K. J., Nave, C. P., & Rasio, F. A. 2001, ApJ, 550, 691 
Joshi, K. J., Rasio, F. A., & Portegies Zwart, S. 2000, ApJ, 540, 
969 

Kalirai, J. S., Ventura, P., Richer, H. B., Fahlman, G. G., Durrell, 

P. R., D'Antona, F., & Marconi, G. 2001, AJ, 122, 3239 
Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 
Mackey, A. D., Wilkinson, M. I., Davies, M. B., & Gilmore, G. F. 

2007, MNRAS, 379, L40 
Marshall, J. R., van Loon, J. T., Matsuura, M., Wood, P. R., 

Zijlstra, A. A., & Whitelock, P. A. 2004, MNRAS, 355, 1348 
Merritt, D., Piatek, S., Portegies Zwart, S., & Hemsendorf, M. 

2004, ApJ, 608, L25 
Noyola, E. & Gebhardt, K. 2006, AJ, 132, 447 
Portegies Zwart, S. F. & McMillan, S. L. W. 2002, ApJ, 576, 899 
Spruit, H. C. 1998, A&A, 333, 603 

Trenti, M. 2006, ArXiv Astrophysics e-prints (astro-ph/0612040) 
Vassiliadis, E. & Wood, P. R. 1993, ApJ, 413, 641 
Weidemann, V. 1977, A&A, 59, 411 

Williams, K. A. 2002, PhD thesis, UNIVERSITY OF 

CALIFORNIA, SANTA CRUZ 
Yungelson, L. R., van den Heuvel, E. P. J., Vink, J. S., Portegies 

Zwart, S. F., & de Koter, A. 2008, A&A, 477, 223 



